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ABSTRACT 

Recent Low Frequency Array (LOFAR) observations at 115 — 175 MHz of a field at 
medium Galactic latitudes (centered at the bright quasar 3C196) have shown striking 
filamentary structures in polarization that extend over more than 4 degrees across the 
sky. In addition, the Planck satellite has released full sky maps of the dust emission in 
polarization at 353 GHz. The LOFAR data resolve Faraday structures along the line 
of sight, whereas the Planck dust polarization maps probe the orientation of the sky 
projected magnetic field component. Hence, no apparent correlation between the two 
is expected. Here we report a surprising, yet clear, correlation between the filamentary 
structures, detected with LOFAR, and the magnetic field orientation, probed by the 
Planck satellite. This finding points to a common, yet unclear, physical origin of 
the two measurements in this specific area in the sky. A number of follow-up multi¬ 
frequency studies are proposed to shed light on this unexpected finding. 

Key words: ISM: general, magnetic fields, structure - radio continuum: general, 
ISM - submillimetre: ISM -- techniques: interferometric, polarimetric 


1 INTRODUCTION 

Non-thermal plasma is one of the main components of the 
interstellar medium (ISM) in our Galaxy (Ferriere 2001). 
It mostly contains relativistic electrons and protons, which 
permeate the ISM and interact with the all-pervasive mag¬ 
netic field of our Galaxy (Haverkorn 2015), resulting in syn¬ 
chrotron emission. At low radio frequencies, like those mea¬ 
sured by the LOFAR telescope (van Haarlem et al. 2013), 
the line-of-sight (LOS) magnetic field component (i3||) and 
ionized gas can be probed by changes in the polarization 
of the Galactic synchrotron radiation through the so-called 
Faraday rotation. 

Recent WSRT and LOFAR observations (Bernardi et 
al. 2009; lacobelli et al. 2013; Jelic et al. 2014; Jelic et al. 
2015) have shown surprisingsly strong polarization at fre- 
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quencies around 150 MHz. This polarization originates in 
the diffuse Galactic background. Due to Faraday rotation 
in the medium between source and observer the frequency 
dependence of the polarization can be very complex. The 
frequency-dependent polarization can be unravelled using 
a technique called rotation measure (RM) synthesis (Brent- 
jens & de Bruyn 2005); this produces images as a function of 
RM or Faraday depth. The observed structures in Faraday 
space can be very complex morphologically. In particular, 
the 3C196 field (Jelic et al. 2015) shows a strikingly straight 
filamentary structure that extends over scales of more than 
4 degrees. This filament also exhibits narrow depolarization 
canals, with a width down to 1 arcmin, that extend along 
the whole structure (Jelic et al. 2015). 

At sub-millimetre wavelengths, in particular at the fre¬ 
quency band of 353 GHz, the Planck satellite measured dust 
polarization with unprecedented sensitivity {Planck Gollab- 
oration XIX 2015; Planck Gollaboration I 2015). Polariza- 
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tion results from the alignment of interstellar dust grains 
with respect to the Galactic magnetic held (Davis & Green- 
stein 1951; Lazarian 2007). Hence, dust polarized emission 
provides powerful insights into the physical properties of the 
dust and the magnetic Held. 

The physical quantities that the two probes measure are 
in general thought to be unrelated. Whereas Faraday rota¬ 
tion at low frequencies is sensitive to localized structures in 
Faraday space, which in the 3C196 Held mostly have small 
Faraday depth, the Planck dust polarization gives the den¬ 
sity weighted mean orientation of the magnetic Held compo¬ 
nent perpendicular to the LOS, Indeed, in a number 

of LOFAR observational windows, e.g., ELAIS-Nl and NCP 
(Jelic et al. 2014), no correlation between the two probes is 
evident. However, in this letter, we report a strong correla¬ 
tion between the two probes in the 3C196 Held. We discuss 
possible further observations that might help understand the 
physical origin of this correlation. 

In § 2 and § 3 of this letter, the LOFAR and Planck 
satellite data towards the 3C196 quasar are presented, re¬ 
spectively. In § 4 we discuss the correlation between the two 
Helds and we conclude with discussion and summary in § 5. 


2 THE LOFAR 3C196 DATA 

The 3C196 Held is one of the primary Helds of the LOFAR- 
Epoch of Reionization key science project (de Bruyn, in 
prep.), centered on the very bright compact quasar 3C196. 
It is located in a radio cold region towards the Galactic an¬ 
ticentre {I = 173°, b = -1-33°). It has been observed on a 
large number of nights with the LOFAR High Band Anten¬ 
nas (HBA). For this paper we use images based on one fully 
processed 8h synthesis, taken in December 2012. A detailed 
description of the data, their processing and the image qual¬ 
ity obtained from them are found in Jelic et al. (2015). For 
completeness, we present here a brief overview of the data. 

The dataset contains all four correlation products be¬ 
tween pairs of orthogonal LOFAR dipoles from which full 
Stokes (IQLfV) images were made in the frequency range 
of 115 — 175 MHz with 183 kHz resolution. The images were 
produced with a frequency independent resolution of about 3 
arcminutes. On these images we apply RM synthesis (Bren- 
tjens & de Bruyn 2005) to study the linearly polarized emis¬ 
sion as a function of Faraday depth. Faraday depth gives the 
amount of Faraday rotation due to the intervening magneto¬ 
ionic medium between the source of the polarized emission 
and the observer, 

^observer 

<I> = 0.81 rad m~^ UeB^^dl. (1) 

^source 

Here rie is the electron number density in cm~^, is mea¬ 
sured in /iG, and dl is the LOS element in pc. The spectral 
resolution and frequency range of the data constrain the res¬ 
olution of our cubes in Faraday depth to = 0.9 rad m~^ 
and the largest Faraday structure that can be detected 
A<F = 1.1 rad (Brentjens & de Bruyn 2005). 

Detected structures of Galactic polarized emission are 
spread over a wide range of Faraday depths (from —3 
to -1-8 rad m~^) and their peak brightness temperature is 
between 5 and 15K (Jelic et al. 2015). Figure 1 shows 
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Figure 1. A composite image of morphological features of the 
3C196 field detected with LOFAR at different Faraday depths and 
the magnetic field lines orientation (gray lines) inferred from the 
Planck dust polarization maps at 353 GHz. A “triangular” feature 
displayed in green (marked with C) is emission at negative Fara¬ 
day depths (—3 to —0.5rad m“^), the filamentary structures at 
Faraday depth of -1-0.5 rad m~^ are given in yellow (marked with 
A and B). The violet shows the prominent diffuse background 
emission arising at Faraday depths from +1.0 to +4.5rad m~^. 
The resolution of the LOFAR image is 3 arcmin. 

the most interesting morphological features of this emis¬ 
sion in a composite image of polarized structures at dif¬ 
ferent Faraday depths. A “triangular” feature showing emis¬ 
sion at negative Faraday depths (—3 to —0.5 rad m~^) is 
shown in green. Relatively straight filamentary structures 
running from South to North at Faraday depths around 
+0.5 rad m“^ are shown in yellow, while the prominent dif¬ 
fuse background emission arising at Faraday depths from 
+ 1.0 to +4.5 rad m~^ is given in violet. The filaments 
are displaced by ss —1.5radm~^ in Faraday depth rela¬ 
tive to the surrounding background emission, suggesting a 
magneto-ionic medium located in front of the bulk of the 
emission (Jelic et al. 2015). Note that Faraday depth is un¬ 
certain at the level of fe 0.1 rad due to the ionospheric 
correction uncertainty. 


3 THE PLANCK MAGNETIC FIELD DATA 

The Planck satellite has measured the sky in 9 frequen¬ 
cies ranging from 30 GHz to 857 GHz. An inspection of 
these maps reveals that the lowest frequencies (30 GHz and 
70 GHz) and the highest frequencies (^ 217 GHz) are dom¬ 
inated by Galactic synchrotron emission and thermal dust 
emission, respectively [Planck Collaboration I 2015). In par¬ 
ticular, the Planck satellite team used the polarized data at 
u = 353 GHz to construct an all-sky map of the Galactic 
magnetic field orientation projected onto the plane of the 
sky {Planck Collaboration I 2015). 

Interstellar dust particles are non-spherical, which gives 
rise to polarized emission refiecting this asphericity. Coher¬ 
ent polarized emission is observed because the elongated 
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dust grains tend to be aligned, along their short axes, with 
respect to the magnetic field. The details of the grain align¬ 
ment process are not completely clear but there is a broad 
agreement that radiative torques are an effective means of 
grain alignment in the diffuse ISM (for more details see 
Draine & Weingartner 1996; Lazarian 2007) are an effec¬ 
tive means of the observed grain alignment in the diffuse 
medium {Planck Collaboration XX 2015). As a result, the 
dust polarization angle is at 90° of the perpendicular com¬ 
ponent of the magnetic field vector averaged along the LOS, 

(B±). 

For this work, we use the line integral convolution (LIC; 
Cabral & Leedom 1993) representation of the magnetic held 
inferred from the 353 GHz polarization maps publicly avail¬ 
able at the Planck Legacy Archive^, smoothed to 15 arcmin 
resolution. LIC is a representation algorithm of a pseudo¬ 
vector held obtained by computing the stream line that 
starts at the centre of each pixel. For this study, we use the 
LIC map corresponding to a maximum integrated length of 
one fourth of the map size per iteration and three iterations 
of the algorithm. These LIC magnetic held lines are shown 
in Fig. 1. The paper results are robust with respect to the 
specihc choice of LIC algorithm parameters. 


4 COMPARISON BETWEEN THE MAPS AND 
ITS INTERPRETATION 

We now inspect the various correlated features between the 
two data sets. There are three regions of interest, marked as 

A, B and C in Fig. 1. The strength of the correlation be¬ 
tween the two probes is discussed separately for each region. 
In general, there is a broad alignment between the South to 
North orientation of the magnetic held lines as measured 
by Planck with the orientation of the LOFAR hlaments de¬ 
tected in Faraday depth. 

In particular, the linearly polarized emission at Faraday 
depth of -1-0.5 rad m~^ at the centre of the image (shown in 
yellow and marked as A) is very strongly aligned with the 
Planck magnetic held orientation. This alignment is not only 
with the feature at this specihc Faraday depth in the LO¬ 
FAR data but also with many other features associated with 
this hlament, such as the depolarization canals described in 
Jelic et al. (2015). This correlation might not be surprising 
since the overall Galactic magnetic held direction, is aligned 
with the South-North direction in this specihc projection. 
However, we note that such a correlation is not seen in the 
other LOFAR observational windows, at comparable Galac¬ 
tic latitudes like the ELAIS-Nl (Jelic et al. 2014) and NGP 
(Jelic et al. in prep.) helds. 

The second correlated feature is shown in the region 

B. Here the bending of the magnetic held line follows that 
of the hlament at the Faraday depth of -1-0.5 rad m”^. The 
correlation of this feature in Fig. 1, though obvious, requires 
a more careful statistical analysis to assess its signihcance. 
In order to quantify this we resort to the RM cubes. We 
use Stokes Q images that contain information about both 
intensity and polarization angle and therefore exhibit larger 
contrast, which helps dehning the features at Faraday depth 

^ http://pla.esac.esa.int 


Stokes Q: $ =+0.5 rad m 



Figure 2. This figure shows in gray scale the Stokes Q image at 
Faraday depth of +0.5 rad (see Jelic et al. 2015, for details). 
The area highlighted with the black contour is the region in which 
the polarization angle is roughly uniform. The red line shows an 
order 10 polynomial fit to the black contour. The Planck mag¬ 
netic field lines are overlaid and indicated with the open yellow 
contours. The magnetic field line shown in green marks the line 
that has a shape similar to the line of uniform polarization angle. 

of +0.5rad m~^ (see Fig. 2). The area deliniated with a 
black contour is the region in which the polarization angle 
is roughly constant (see Jelic et al. 2015, for more details). 
The Planck magnetic field lines are overlaid and indicated 
with the yellow contours. The filled line in green marks the 
Planck magnetic field line that has a shape similar to the 
LOFAR feature. 

In order to define these features both are fitted with or¬ 
der 10 polynomials, resulting in the red line for the LOFAR 
data shown in the hgure^. The angle, 6, of the tangent at 
each point along the curve is calculated for both LOFAR 
(red) and Planck (green) lines and compared using a nor¬ 
malized x-squared statistic. For two lines, or part of them, 
to be similar the angles defining their direction should be 
close over a long stretch of their path. Figure 3 shows the 
square of the difference between these angles, which we call 
9^ and 9^ for LOFAR and Planck, respectively. The red 
curve in Fig. 3 corresponds to the squared angle difference 
between the most similar lines in Fig. 2 (red and green lines), 
whereas the others correspond to the squared angle differ¬ 
ence between the red line and hve randomly chosen Planck 
magnetic held lines, also from Fig. 2. Clearly, the red curve 
shows the smallest difference for the longest stretch of length 
along the two curves. The two vertical dashed lines indi¬ 
cate the range over which the squared angle difference is the 
smallest (< 10~^ for about 75 arcmin). 

^ The choice to fit order 10 polynomial is somewhat arbitrary. It 
serves only to fit a smooth line to the data that allows a clear 
determination of the orientation of the magnetic and constant 
polarization angle lines. The only requirement is that the order 
of the polynomial should neither be too small nor too large, so as 
not to underfit or overfit the data, respectuvely. 
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Figure 3. This figure shows the square of the difference between 
the angle of the tangent at each point along the LOFAR and 
Planck curves, 0^ and 9^, respectively. The red curve in the 
figure shows the squared difference between the angles of the most 
similar lines in Figure 2 (red and green lines), whereas the rest 
shows the angle difference squared between the yellow line and five 
randomly chosen Planck magnetic field lines, also from Figure 2. 
The two vertical dashed lines mark the range over which the red 
line is the smallest; namely < 10~^ for about 75 arcmin). 

Therefore, we define our statistic to be 

= ^ T. idt-eff, ( 2 ) 

where Npij; is the number of adjacent pixels aronnd the 
pixel with the minimum angle difference; is the angle 
of the LOFAR curve tangent at pixel i and Of is that of 
Planck. Npix is chosen to be 150 pixels which corresponds 
to 75 arcmin, defined by the two dashed lines in Figure 3. 
We note that since the two curves have in general differ¬ 
ent lengths, the comparison is carried out within the over¬ 
lap region. In case this overlap covers less than 150 pixels, 
we choose Npi^ = min{150, where Noveriap is the 

number of pixels in the overlapping region. 

To address how unusually similar these two curves are, 
all the Planck magnetic field lines apparent in Fig. 2 (about 
40) are compared to the LOFAR structure of region B us¬ 
ing the statistic defined in Equation 2. To obtain more ran¬ 
dom realizations of the magnetic field lines that retain their 
statistical and physical properties we randomly shift the 
Planck field lines along the Declination, and if needed, we 
wrap them across the Declination boundaries. This results 
In about 1000 comparisons. The probability density func¬ 
tion of this statistic is shown in Figure 4 with logarithmic 
bins. The left side of the histogram is htted with a Gaussian 
(the solid smooth line). As expected the PDF is asymmetric 
around its peak, however, we note that the existence of a 
clear sharp peak is due to the special nature of the 3C196 
held in which most features are along the South-North di¬ 
rection. The red vertical line in this hgure shows the result 
of the statistic for the most similar lines in Fig. 2 (red 
and green lines). This clearly shows that the similarity be¬ 
tween these two lines is statistically signihcant and unlikely 





Figure 4. This figure shows the probability density function of 
statistic defined in Eq. 2 with logarithmic bins. The left side of 
the histogram is fitted with a Gaussian (the solid line). This is 
in practice a log-normal distribution. The Gaussian is centered 
at -1.38 and its width is ctq = 0.5. The red vertical line shows 
the value of the statistic comparing the most similar lines in 
Figure 2 (red and green line). 

due to statistical coincidence. Notice that this statistic does 
not take into account the location of each curve, but rather 
compares only their shape. Therefore, the fact that these 
two curves are unusually similar and are co-located on the 
sky makes the statistical significance of this correlation even 
larger. We note that this result was confirmed by similar 
analysis based directly on the Planck Q and U maps, inde¬ 
pendently of the Lie algorithm. 

The third region that shows correlated features is re¬ 
gion C in Figure 1. Around this region the Planck magnetic 
field lines broadly follow the bend in the red line as indi¬ 
cated by the arrow in Figure 2. Moreover, the magnetic field 
lines from Planck seem to be less coherent within this re¬ 
gion. We have deconvolved the magnetic field into rotational 
(curl) and potential (divergence) components, assuming the 
third direction to be the LOS direction. This analysis shows 
that the magnetic field within region C have a dominant 
rotational component. This component can not be due to 
turbulence simply because region C has a coherent Faraday 
displacement. The last two statements are hard to quantify 
statistically, hence they should be viewed merely as poten¬ 
tially interesting observations. 


5 DISCUSSION AND CONCLUSIONS 

In the previous section we have discussed the morphology 
and significance of the correlated features between the LO¬ 
FAR and Planck data sets in the 3C196 field. The LO¬ 
FAR polarization data in Faraday depth is proportional to 
neBji Al (Eq. 1), whereas Planck measures, through the po¬ 
larized dust emission, (By). Since dust and electron distri¬ 
butions are probably unrelated and the magnetic field com¬ 
ponents that these probes measure are different, one does 


MNRAS 000, 000-000 (0000) 










































Galactic interstellar filaments as probed by LOFAR and Planck 5 


not generally expect the two data sets to be correlated. In¬ 
deed, in a number of common observational fields, we con¬ 
firm that such a correlation is not evident. 

The correlation features reported in this paper point 
to a common underlying physical structure that influences 
both the LOFAR and Planck data sets. One can roughly 
divide the features into global and local features. The global 
feature is the broad correlation of the dust inferred orienta¬ 
tion of the magnetic field and the LOFAR filaments, clearly 
apparent in the images. This is especially striking when one 
inspects features A and B in Figure 1. Jelic et al. (2015) 
make the point that since in the 3C196 field the magnetic 
field that follows the spiral arm of our Galaxy is almost per¬ 
pendicular to the LOS, i3|| is small relative to {B±). This 
conclusion is supported by the estimated {B||) = 0.3+0.1 pG 
in the direction of the pulsar J081558+461155, located about 
2° south of 3C196, using the RM to dispersion measure 
ratio (see Jelic et al. 2015, and references therein). A low 
J3|| suggests that the origin of the correlation resides in the 
alignment of the electron density distribution with {Bj_). 
Such alignment has been observed for the filamentary struc¬ 
tures traced by dust emission (Glark, Peek, & Putman 2014; 
Planck Gollaboration Int. XXXII 2015). In other words, the 
broad alignment between the two images might be driven by 
the global magnetic field of the Galaxy, which is typically 
~ 5 fj.G (Haverkorn 2015). This interpretation, although the 
most straightforward for us, does not preclude other expla¬ 
nations that rely on variations in the structure of i3||. 

The correlations implied by the local features are appar¬ 
ent in regions B and C. The agreement between the bend¬ 
ing of the lines in both regions, especially in B, argues that 
the two effects come from the same physical location. The 
physical properties of the ISM in this region are uncertain. 
However, some constraints can be derived from the LOFAR 
observations (Jelic et al. 2015). The lack of observed emis¬ 
sion in the LOFAR total intensity, associated with the fila¬ 
mentary structures detected in polarization, gives an upper 
limit to their thermal free-free emission, Ta < 0.2 K. If we 
adopt Te = 8000 K for the electron temperature (Reynolds 
1991), which is a typical value for the warm ionized ISM, 
and Al = 1 pc for the LOS thickness (Jelic et al. 2015), 
we find Ue ^ lcm~^. The thickness of the structures in 
Faraday depth. A"!* ~ lradm~^, constrains i3|| > 1.2/rG, 
assuming again Al = 1 pc. Motivated by the agreement in 
the bending of region B, we consider equipartition between 
magnetic energy and thermal energy. Using these constraints 
on Us, we find that the strength of the total magnetic field 
has to be Btotai ^ 6.5 ^G. Note that in reality these values 
might be different as equipartition might not hold between 
the magentic energy and thermal energy but rather between 
the magnetic and kinetic energy associated with bulk mo¬ 
tions. We also note that the true values of J3|| and Btotai are 
probably close to the limits derived above, because obvi¬ 
ously i3|| Btotai, moreover the detection of polarized dust 
emission by Planck requires a significant {Bj_). This is also 
consistent with the analysis of the LOFAR data in this field 
as noted by Jelic et al. (2015). 

To resolve the ambiguity of the physical parameters and 
shed light on the physical picture of the ISM, additional 
observational probes are needed. Here we mention a few. 

• Dust particles are usually embedded within neutral hy¬ 


drogen (Hi) gas. This is especially evident through Hl/dust 
correlation at high Galactic latitudes, where the sky is 
largely transparent to Hi emission and infrared emission 
from dust. Therefore, Hi data can provide kinematic infor¬ 
mation about the features studied in this paper. 

• Observed starlight polarization, caused by selective ab¬ 
sorption by magnetically aligned interstellar dust grains 
along the LOS, traces the orientation of (Bx) (Hiltner 1956; 
Serkowski, Mathewson, & Ford 1975; Pavel 2014). By re¬ 
peating this measurement for stars at different distances, 
one can map the the 3D distribution of the orientation of 
(B ±) and provide complementary information to Planck's 
dust polarization observations. However, this type of po- 
larimetry is sensitive only to the orientation of the interstel¬ 
lar magnetic field. It is worth noting that a 3D mapping of 
the Galactic magnetic field might be crucial for the removal 
of dust polarisation foreground in order to measure the pri¬ 
mordial GMB B-mode polarisation signal (see e.g., Tassis & 
Pavlidou 2015). 

• To determine the strength of the interstellar magnetic 
field, one can attempt to measure the Zeeman splitting of 
the Hi 21-cm line or other radio lines. Such data measure the 
Bjl component of interstellar magnetic field and complement 
observations of the Faraday depth structure. 

• Future limits on the fine structure in total intensity of 
the diffuse emission will allows us to constrain |Bxp. 

Finally, the correlation between these two probes calls 
for exploration of a similar correlation in other fields. Specif¬ 
ically, we propose to explore regions in which the Planck 
derived magnetic field is coherent and ordered over a large 
region of the sky and follow such regions with LOFAR ob¬ 
servations. This comparison might shed more light on the 
physical origin of this correlation. 
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